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' Abstract In the last decade, the photospheric solar metallicity as determined 

$H . from spectroscopy experienced a remarkable downward revision. Part of this 

^ ' effect can be attributed to an improvement of atomic data and the inclusion 

of NLTE computations, but also the use of hydrodynamical model atmospheres 
l/^ ', seemed to play a role. This "decrease" with time of the metallicity of the solar 

photosphere increased the disagreement with the results from helioseismology. 
■ With a C05B0LD 3D model of the solar atmosphere, the CIFIST team at the 

, Paris Observatory re-determined the photospheric solar abundances of several 

(-H ' elements, among them C, N, and O. The spectroscopic abundances are obtained 

CIh' by fitting the equivalent width and/or the profile of observed spectral lines with 

^ , synthetic spectra computed from the 3D model atmosphere. We conclude that 

^ ' the effects of granular fluctuations depend on the characteristics of the individual 

^ ■ lines, but are found to be relevant only in a few particular cases. 3D effects are not 

reponsible for the systematic lowering of the solar abundances in recent years. 
The solar metallicity resulting from this analysis is Z = 0.0153, Z/X = 0.0209. 
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1. Introduction 

Since the pioneering work of IRusselll (|1929p . the study of the chemical composi- 
tion of the solar photosphere has been an important topic in astronomy. Russell, 
measuring the strength of the absorption lines in the observed solar spectrum, 
determined the abundance of 56 elements and six molecules. After this work 
many other astronomers analysed the solar spectrum in order to deduce the 
detailed pattern of photospheric abundances. The accurate determination of the 
solar chemical abundances is a major topic, because: 

1. the knowledge of the elemental abundances in the present solar photosphere 
is the basis to infer the chemical composition of the initial Sun, and allows one 
to reconstruct the past and future evolution of the Sun, including its physical 
and chemical internal structure; 

2. the comparison of the chemical abundances in the solar photosphere and in 
meteoritic samples provides important information about the formation and 
chemical evolution of the solar system; 

3. the construction of detailed models of the solar atmosphere requires the 
knowledge of its chemical composition; 

4. the knowledge of the solar photospheric chemical abundance allows the em- 
pirical determination of the oscillator strength of any spectral line of this 
element observable in the solar spectrum; 

5. the solar abundances serve as a reference for the chemical analysis of other 
stars in the Galaxy, of the interstellar medium, and of the stellar populations 
of external galaxies. 

In our analysis of the solar photosphere we are mainly interested in the first and 
third point. 

For our work on the solar abundance determinations, we rely on a CO^BOLD 
3D model of the solar photosphere. We were interested to understand if the 
presence of horizontal fluctuations (the solar granulation) has a systematic effect 
on the abundances derived from the 3D model with respect to what is obtained 
by ID models which ignore granulation. 

We find that the effects of granular fluctuations can lead to both positive and 
negative abundance corrections, depending on the properties of the individual 
spectral line under consideration. However, the granulation effect is relevant 
only in a few particular cases. 3D effects are not responsible for the systematic 
lowering of the solar abundances in recent years. 

In the following sections we summarise the investigations of some elements we 
have already published, and present a new analysis of some other elements (Li, K, 
Fe, and Os) . For a complete overview of the solar abundance determinations we 
recommend the review bv lLodders et aLKpOO^ as a very detailed and complete 
work. 

2. Model Atmospheres 

Our abundance analyses of the solar photosphere are based on a 
time-dependent, 3D, hydrodynamical model atmosphere computed with 
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the CO^BOLD code |http://w ww.a5tro. uu.se/ -bf/co5bold-20020216/cobold.htm 



dFreytag, Steffen, andDorch| [20021 12003) • The 3D model atmosphere we use 
has a box size of 5.6 x 5.6 x 2.27 Mm'^, a resolution of 140 x 140 x 150 grid 
points, and spans a range in optical depth of about —6.7 < logrRoss < 5.5 
(from —1.4 Mm below to +0.9 Mm above tross = !)• The selected 19 snapshots 
wc use for the spectrum-synthesis calculations arc equidistantly spaced in time, 
sufficiently separated in time to show little correlation, and cover 1.2 hours of 
solar time. For comparison, the following ID plane-parallel model atmospheres 
were considered: 

1. the (3D) model obtained by horizontally averaging each 3D snapshot over 
surfaces of equal (Rosseland) optical depth; 

2. the IDlhd model, a ID hydrostatic mixing-length model which employs the 
same micro-physics and radiative transfer scheme as CO^BOLD; 

3. the semi-empirical Holweger-Miiller solar model ( |Holweger| 119671 
[Holweger and Miiller|[Wl hereafter HM); 

4. the ATLAS-9 model computed by F. Castelh 



^http://w wwuser. oats, inaf.it/castelli/sun/ap00t5777g44377klasp.dat with 
the solar abundances of |Asplund, Grevesse, and Sauval| (|2005p . 

The (3D) and IDlhd models have been introduced as reference 
ID models, because they share the microphysics with the 3D model 
( [Caffau and Ludwig, 2007 ). For more details about the ID and 3D models used 



in our analysis see ICaffau et al.\ (|2008ap or the other papers of the collaboration 
on the determination of the solar photosphcric abundances. 



3. 3D Corrections 

We investigated the effects of 3D convection on solar abundances in 
[Caffau, Ludwig, and Steffen] (j2009p . To define 3D corrections we selected two 
reference models: (3D) and IDlhd- Both share the microphysics of the 3D- 
CO^BOLD model to ensure differential comparability. We distinguish two 
different 3D corrections, defined as: A^^\^^ic) = A{Y)3u — ^(^){3D> and 
mic 5 ttMLx) = ^(^)3D — ^(^)iDlhd: where Y represents any chemical ele- 
ment, and A{Y) = log (ny /jih) + 12. The first correction isolates the granulation 
effects, i.e. the infiuence of the horizontal fluctuations around the mean strati- 
fication, while the second one measures the total 3D effect, accounting for both 
horizontal fluctuations and for the different influence of 3D hydrodynamical and 
ID mixing-length convection, respectively, on the resulting mean temperature 
structure. Obviously, both 3D corrections depend on the choice of the microtur- 
bulence parameter, [^mic]i used with the ID models, while the second one is also a 
function of the mixing length parameter, [aMLx], adopted for the IDlhd model. 
The idea is that these 3D corrections can be derived from the simulations with 
much better accuracy than the detailed thermal structure of the solar atmosphere 
in absolute terms. While the latter depends sensitively on e.g., the frequency 
binning adopted for the radiation hydrodynamics, the 3D corrections are less 
sensitive to such details, as they are defined in a strictly differential way. Hence, 
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the 3D corrections are not merely evaluated to measure the abundance difference 
between a given pair of 3D and ID models, but they may rather be utilised to 
improve the abundance determinations from standard ID models. Abundances 
derived from the semi-empirical HM model should be corrected by adding only 
the "granulation correction" A^^'' (^mic), while abundances based on standard ID 
mixing- length models like ATLAS or MARCS, must be corrected by adding the 
"total 3D correction" A'^^(^inic, aMLx)- 

The 3D corrections are well defined as long as the considered lines are weak 
such that their equivalent widths are independent of the assumed microtur- 
bulence. The situation becomes more problematic once the lines are partly 
saturated, because then the corrections depend on ^mic, and the choice of the mi- 
croturbulence parameter is critical. In principle, the microturbulence parameter 
derived empirically from analysing observed solar spectra with a standard ID 
model, ^mic(Sun-lX). should be identical with the theoretical microturbulence 
?mic(Hydro-lDLHD)j obtained from analysing the synthetic spectrum of the 3D 
hydrodynamical model atmosphere with the IDlhd model. However, there are 
indications that ^mic(Hydro-lDLHD) is systematically smaller than ^inic(Sun-lX) 
( Steffen et al, 2009] ). As a consequence, the 3D corrections should be computed 
with ^mic (Hydro- IDlhd)- Wherever possible, abundances from weak lines are to 
be preferred. 

For some of the weak lines that we investigated {e.g. the Ci line at 538 nm), 
we compared the 3D correction A^^^ with the results of [Steffen and Holweger| 
(|2002p . which are based on a 2D hydrodynamical simulation, and find a close 
agreement. Summarising our analysis, we can say that the total 3D correction 
A*^^) is positive for the majority of the relevant spectral lines, except for a few 
weak, high-excitation lines (mostly Nitrogen) for which both A*^^^ and A^^^ are 
negative. However, the 3D corrections are small in general. 

The abundances presented in this work have been derived directly from the 
3D model, in line with our previous publications. We trust in the temperature 
structure of our 3D-C0^B0LD model, because it is able to reproduce the centre- 
to-limb variation of the continuum intensity even somewhat better than the HM 
model (Ludwig et al.., 2009). However, we plan to compute abundances also with 
the alternative approach, i.e. from the HM model and correcting by A(^)(^niic), 
in future investigations. While it is not entirely clear which of the two approaches 
gives the more reliable results, we propose the abundance differences between 
the two methods to be considered as a measure of the systematic uncertainty of 
the abundance analysis. 



4. Observational Data 

Our analysis is based on mainly four high resolution, high signal-to-noise ratio 
[S/N] spectra, two for disc centre, and two for the integrated disc. For most 
elements we use more than one solar spectrum, because we realised that the 
abundance derived from different spectra do not always agree within one a. The 
observed spectra we considered are: 
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1. the integrated disc spectrum based on fifty solar FTS scans taken by J. Brault 
and L. Testerman at Kitt Peak between 1981 and 1984 ( |Kurucz, 2005D ; 

2. the two absolutely calibrated FTS spectra obtained at Kitt Peak in the 1980s, 
covering the range 330 nm to 1250 nm for the integrated disc and disc centre 
dNeckel and Labs, 198"4| [Neckel, 1999] ); 

3. the disc centre intensity spectrum in the range 300 nm to 1000 nm ob- 
served from the Jungfraujoch ( [Delbouille, Roland, and Neven, 1973] ), and 
in the range from 1000 nm to 5400 nm observed from Kitt Peak 
Pelbouille el at, 198T| ). 

5. Chemical Abundances 
5.1. Lithium 

Lithium is widely studied in metal-poor stars. The Lii resonance doublet at 
670.7 nm is also observable in the solar spectrum and in the spectra of F-K 
main sequence stars. In metal-poor stars, this region is very clean, and the 
Lithium feature is not blended. This is not the case in solar-metallicity stars, 
where atomic and molecular lines contaminate the region. Good atomic data for 
these blending lines arc not readily available in databases, but there are a few 
published line lists that have been used for the abundance analysis of the solar 
photosphere. For our own analysis of Li in the solar photosphere, we took into 
account both granulation and 3D-NLTE0 effects in computing the contribution 
of the Li doublet, while 3D-LTE line formation calculations were performed for 
the other blend components. 

Considering the four solar atlases described above to investigate lithium, we 
realise that the two observed intensity spectra are significantly different (see 
Figured]). We have no explanation for this disagreement. The easiest way to 
explain it would be to invoke telluric absorption, but this spectral range has 
been carefully scrutinised in the context of the Li isotopic ratio in metal-poor 
stars, and it seems unlikely that a telluric absorption could have gone unnoticed. 
We decided to discard the Delbouille disc centre atlas, and to work only with the 
other three atlases for the Li abundance determination. The reason for discarding 
the Delbouille disc centre atlas is that we cannot reproduce the profile with 
synthetic spectra, while we can in the case of the other three solar atlases, and 
obtain Li abundances that agree closely. 

The two lists of blending lines that can reproduce the solar spectrum in this 
region reasonably well are those from lGhezzi et aL\l\2009\\ and from |Reddy et aL\ 
(|2002p . To obtain the Li abundance, we fit the observed line profile, interpolating 
in a grid of 3D synthetic spectra, in which NLTE effects are taken into account 
for Li. We obtain A(Li)3D-NLTE = 1-03 ± 0.03, fitting the three observed line 
profiles using the IGhezzi et aL\ (|2009p line list. Fortunately, the solar Lithium 
abundance is not very sensitive to the choice between the two line lists: when 
using the |Reddy et flL|(|2002p list, we obtain A(Li)3D-NLTE = 1-01. 



NLTE stands for Non Local Thermodynamic Equilibrium 
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5.2. Carbon 



For a selection of 45 lines, we measured the equivalent widths, taking into account 
the blending components, when present, and we obtained A{C)= 8.50 ± 0.06 
Caffau et al., 20^. NLTE co rrections were computed with the Kiel code 



Steenbock and Holweger, 19841, using the (3D) model as input. In our sample 



of Ci lines there is one forbidden line at 872.7 nm, formed in LTE conditions. 
Its abundance of A(C)=8.48 ± 0.02 is very close to the value we obtain from 
the complete sample. A considerable fraction of our selected lines are strong, in 
the saturated part of the curve of growth. We do not see a trend of the Carbon 
abundance with respect to the equivalent widths of the lines, and the A(C) is 
unchanged if we consider only lines whose EW is smaller than 12 pm. When we 
consider lines weaker than 8 pm of equivalent width, A(C) is lower of 0.02 dex. 



5.3. Nitrogen 



Our analysis relies on two sets of equivalent width measurements 
( [Biemont et al., 1990| |Grevesse et al., 1990[ ). The Nitrogen abundance is ob- 
tained by interpolating in synthetic curves-of-growth from 3D and ID models. To 
account for NLTE effects, we applied the corrections computed in using the (3D) 
as the ID background stratification. The result is A(N)= 7.86±0.12, as reported 
in detail bv lCaffau efaDipgng)) and lMaiorca et aZ.I([MDg| . [Asplund etld] ([200^)) 
criticise our line selection, but unfortunately provide no information on their own 
analysis, except their line-to-line scatter of 0.04 dex. We note that the difference 
between the abundance from their atomic and molecular NH Aw = 1 lines is 
about 0.1 dex. 



5.4. Oxygen 



Our photospheric Oxygen abundance of 8.76 ± 0.07 is derived from ten atomic 
Oxygen lines. The NLTE corrections arc computed from the (3D) model as 
input to the Kiel code, and applied to the 3D-LTE abundances. For details 
of the analysis see ICaffau et al.\ (|2008ap and [Caffau and Ludwig] (j2008p . In our 
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analysis we find a disagreement between the two forbidden [01] lines at 630 nm 
and 636 nm of more than 0.1 dex that we are not able to explain, unless we 
allow the contribution of the Nickel blend of the line at 630 nm to be smaller 
than expected by about a factor of two (see also [Ayresl (2008)). Interestingly, 
ILambertl p978p already considered the presence of the Ni i line, but concluded 
that such contribution must be small to have an agreement between the two 
[01] lines. This disagreement is not present in |Asplund et al^ (|2004p . while in 
[Asplund et al.\ (|2009p a disagreement is recognised, but no solution is offered. 
When we use the equivalent widths given by [Asplund et a/](|2004p together with 
our 3D model, we find indeed an agreement between the two lines, but with a 
"high" A{0) value of 8.79 and 8.76, respectively, while |Asplund et a/.|(PM)) give 
A{0) = 8.69 and 8.67 for the two forbidden lines. This difference is probably 
related to differences in the temperature structure of the two independent 3D 
model atmospheres. 



5.5. Phosphorus 



We considered five infrared P i lines of Multiplet 1 for the photospheric abun- 
dance determination ( [Caffau et ai, 2007bp . Our photospheric Phosphorus abun- 
dance of A{F)= 5.46 ± 0.04 is in perfect agreement with the meteoritic value 
dLodders et ai, 2009p . The Phosphorus abundance in [Asplund et al.\ (|2009p 
is 0.06 dex larger with respect to the value of A{P)= 5.35 ± 0.04 given in 
[Asplund, Grevesse, and Sauval (|2005p . It is unclear whether this upward revision 
is related to their new 3D model or to a change in the line list, in the log gf 
values, and/or the adopted equivalent widths. 



5.6. Sulphur 

In lLodders et aLI(|2009p the photospheric Sulphur abundance is given as A(S)= 
7.14±0.01 and the meteoritic one as A(S)= 7.17±0.02. We have studied several 
S I lines in the solar spectrum. From the weak forbidden line at 1082 nm we obtain 
A{S)= 7.15 ± (0.01),t^t ± (0.05),y, dCatfau and Ludwig, 20071 ). [Asplund et aL\ 
(j2009p criticise both the measured equivalent widths as being too small, and the 
log gf we (and also [Ryde| (|2006p ) use as being obsolete, claiming that updated 
values should be used. As no further details are given, it is impossible to judge 
whether their equivalent widths and log gf values are better. We considered 
also the permitted lines of Multiplet 3, Multiplet 6, and Multiplet 8, discarding 
the strong lines of Multiplet 1 at 920 nm, because they are blended with telluric 
absorption ( [Caffau et al., 200"7a| . Both lines of Multiplet 6 and Multiplet 8 are 
weak and close to LTE. The abundance we find is ^(S)= 7.14 from the line 
of Muhiplet8, and ^(S)= 7.11 from the two fines of Multiplet 6. The lines of 
Multiplet 3 at 1045 nm are affected by departures from LTE. If we take this 
effect into account, using the NLTE correction of lTakeda et aLI(|2005p . we obtain 
A(S)= 7.30. This value becomes ^(S)= 7.28 when using the NLTE corrections 
computed by S. Andrievsky and S. Korotin with the Sulphur model atom de- 
scribed in IKorotinI ([^005)) and IKorotiiil (1200^)) . The simple average over all of 
the multiplets would give yl(S)= 7.17 ± 0.07, where a is the root-mcan-square 
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deviation. Giving twice the weight to the [SI] lines and the hnes of Multiplet 8 
because they are unaffected by NLTE, we obtain A(S)= 7.16 ± 0.05, which we 
recommend as the solar photospheric value. 

5.7. Potassium 

In the solar spectrum, Potassium is observed through lines of the neutral species, 
K I. As for other alkalines, the strongest lines are those of the resonance doublet 
of Ki at 766.4 nm and 769.8 nm, the only observable lines in metal poor stars. 
In fact, the other Ki lines are much weaker. Unfortunately, the strongest Ki 
line at 766.4 nm is heavily blended in the solar spectrum by strong telluric O2 
absorption. 

Ide La Reza and Miillerl (I1975|) performed the first detailed study of the K i 
769.8 nm doublet component in the solar photosphere, analysing its centre-to- 
limb variation and computing the K abundance in NLTE. 

Among the few excited Ki lines, the one at 404.4 nm is often used in 
old solar abundance analyses; but this line is significantly blended and the / 
value is uncertain so that it is excluded by Ide La Reza and Miillerl (|1975|) . Also 
ILambert and Luck! p978p suggest that a wrong / value could be the reason for 
the exceedingly low abundance derived from this line. 

ITakeda et al.\ Ijl996p obtain a good fit of the solar flux spectrum with their 
synthetic profile of the line at 769.8 nm with a multi-parameter fitting method. 
They also studied the dependence of the K abundance upon the atomic parame- 
ters, the microturbulence, and the adopted model atmosphere (see their Table 4). 
The NLTE correction for this line is estimated by them to be about —0.4 dex. 

|Zhang et al.\ (|2006p consider seven lines: the line at 693.8 nm, the resonance 
doublet, plus four lines in common with ILambert and Luckl (|1978p . These au- 
thors derive the Potassium abundance by line profile fitting of the integrated disc 
solar spectrum with a NLTE synthetic profile. The NLTE corrections provided 
in this paper are applied to our LTE analysis (see below). ITakeda et al. \ (|1996p 
select eight weak K lines, adding to the list of | Zhang et al.\ (|2006p the lines 
at 404.7, 533.9, and 583.1 nm, while they igno re the 1243.2 nm line, w hich is 



included by [Zhang et al.\ (|2006p . According to llvanova and Shimanskiil (|2000p 
NLTE corrections are small for the 693.8 nm line, important for the 1252.2, 
1243.2, and 1176. 9nm lines, and very important for the 769. 8nm line. 

Our present analysis is based on the six K I lines listed in Table[TJ By means of 
the IRAF task splot | http://iraf.noao.edu/ , we have measured the equivalent 
widths in the observed solar atlases. Interpolating in the curve-of-growth, we 
determine the Potassium abundance in the solar photosphere from both 3D and 
ID model atmospheres. NLTE corrections were taken from Zhang et al.\ (j2006p 



where available. The NLTE K abundance derived from the IDlhd model (with 
^mic = l.Okms"^), A{K)= 5.06 ± 0.04, is in good agreement with the value 
of ^(K)=5.12 ± 0.03 obtained by [Zhang etal] ((2006P by using the Kitt Peak 



integrated disc atlas, and excluding the uncertain 766.4 nm line. The Potas- 
sium abundance given by ILodders et t^IUU^ is A(K)= 5.12 ± 0.03, while 
[Asplund et Td] iMM\ recommend A(K)= 5.03 ± 0.09. 



Our adopted value is A(K)3d-nlte = 5.11 ± 0.09, which is obtained by 
applying the NLTE corrections of [Zhang et al.\ (I2006P to the results from the 
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Table 1 . Atomic data and derived abundances for our selection of K I lines 



A 


X 


log 9f 






A(K) 










3D-LTE 


(3D)-LTE 


IDlhd-LTE 


3D-NLTE 


nm 


eV 












580.1749 


1.617 


2.20 


5.263 


5.279 


5.225 


5.203 


693.8763 


1.617 


4.00 


4.998 


5.017 


4.962 


4.929 


769.8974 


0.000 


170.00 


5.434 


5.458 


5.353 


5.144 


1176.9689 


0.000 


50.00 


5.329 


5.359 


5.288 


5.199 


1243.2273 


1.617 


56.00 


5.279 


5.308 


5.237 


5.129 


1252.2134 


1.610 


85.00 


5.275 


5.304 


5.222 


5.055 


(A(K)) 






5.263 


5.296 


5.223 


5.110 



3D-LTE analysis for the integrated disc spectrum. The 3D-LTE resuhs from disc 
centre and integrated disc agree very closely, the latter being 0.03 dex higher. 
This tiny difference is probably due to different NLTE corrections. 

5.8. Iron 



Iron is a complex atom with a very rich spectrum of atomic lines, especially 
from the neutral and singly ionised stage. Many scientists devoted time to the 
determination of the Iron abundance in the solar photosphere, both from Fei 
(which accounts for most of the lines in the solar spectrum) and from Fe ii (which 
is the dominant ionisation stage of Iron in the photosphere) . 

The Iron abundance was continuously decreasing from the value given by 
IRusselll (|1929|) until the 1960s. The correction of the transition probabilities 
produced an increasing solar abundance by about one order of magnitude until 
the 1990s. Since then, A{Fe) has experienced a smooth decrease (see Figure 1 
in lGrevesse and Sauvall (|1999p ). We recall some important developments in this 
field: Corliss and his group investigated the Iron abundance, with their transition 
probability ( [Corliss and Warner, 1966] [Corliss and Tech, 1968[ ); Bridges and 
coworkers applied their new gf-scale to the abundance determination 
( [Bridges and Wiese, 1970| [Bridges and Kornblith, 1974[ ). Between 1980- 
1990, there was a long debate between Holweger and the Blackwell group, 
advocating a low {A{Fc)— 7.50) and a high {A{Fc)— 7.63) Iron abundance, 
respectively ([Biemont et ai, 1991] 



Holweger, Kock, and Bard, 1995 



Blackwell, Lynas-Gray, and Smith, 1995 



Blackwell, Smith, and Lynas-Gray, 1995 



Kostik, Shchukina, and Rutten, 1996] ) . Kostik, Shchukina, and Rutten (I1996P 



pointed out the necessity to perform a 3D-NLTE analysis, which was impossible 
at the time for the lack of computer power and atomic-physics data. 

To investigate the Iron abundance in the solar photosphere, we selected 15 
Fell lines and measured their equivalent widths with the IRAF task splot (see 
Table[3|). We presume that transitions from the dominant ionisation stage of 



Iron, Fell, are close to LTE condition (Gehren et 



for Fei, for which NLTE-effects are important, at the level of 



, 2001 )■ This is not the case 
0.1 dex, as 
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Table 2 . Atomic data and equivalent widths for our 
selection of 15 Fell lines 





X 


log 


nf 


EW [pm] 




eV 


H 


MB 


I 


F 


457.6340 


2.84 


-2.94 


-2.95 


6.90 


6.90 


462.0521 


2.83 


-3.21 


-3.21 


5.60 


5.50 


465.6981 


2.89 


-3.59 


-3.60 


3.60 


3.50 


523.4625 


3.22 


-2.23 


-2.18 


8.80 


8.60 


526.4812 


3.23 


-3.25 


-3.13 


4.70 


4.70 


541.4073 


3.22 


-3.50 


-3.58 


2.90 


2.80 


552.5125 


3.27 


-3.95 


-3.97 


1.35 


1.25 


562.7497 


3.39 


-4.10 


-4.10 


0.82 


0.84 


643.2680 


2.89 


-3.50 


-3.57 


4.40 


4.35 


651.6080 


2.89 


-3.38 


-3.31 


5.80 


5.80 


722.2394 


3.89 


-3.36 


-3.26 


1.95 


1.92 


722.4487 


3.89 


-3.28 


-3.20 


2.30 


2.03 


744.9335 


3.89 


-3.09 


-3.27 


1.75 


1.80 


751.5832 


3.90 


-3.44 


-3.39 


1.60 


1.50 


771.1724 


3.90 


-2.47 


-2.50 


5.30 


5.12 



H: IHannaford et 0^.111199211 
MB: |Melendez and Barbuy| ll2009ll 



shown by IGehren et al.\ (|20011 1200ip using hydrostatic ID model atmospheres, 
and by |Shchukina and Trujiho Bueno| (|200ip using a single snapshot of a 3D 
hydrodynamical simulation. 

For the disc centre and integrated disc spectra, the 3D abundance we find 
is A{Fe)= 7.525 ± 0.057 and A{Fe)= 7.512 ± 0.062, respectively, when using 
the log gf from Melendez and Barbuy| (j2009|) . When wc use the log gf from 
IHannaford et al.] (|1992p . the line-to-line scatter increases and the abundances 
become A{Fe)= 7.530 ± 0.110 and yl(Fe)== 7.516 ± 0.109, respectively If wc re- 
move the 744.9 nm line, which shows an exceptionally low A(Fe) value, we slightly 
reduce the cr: A{Fc)= 7.550±0.082 and A{Fc)= 7.532±0.093, respectively. Thir- 
teen of our selected Fe ii lines are in common with those for which [5chnabel et aL\ 
(|1999|) have measured oscillator strengths. Using their set of log gf values we 
obtain, for these 13 lines, A(Fe)3D = 7.50 ± 0.11 and A(Fe)3D = 7.49 ± 0.11 
for the for the disc centre and integrated disc spectra, respectively. We con- 
clude that the different sets of oscillator strengths provide highly consistent 
Iron abundances. However, the line-to-line scatter is considerably smaller when 
adopting the [Melendez and Barbuy| (I2009P log gf values. Therefore, considering 
both disc centre and integrated disc spectra, we recommend A{Fe)= 7.52 ± 0.06 
for the photosphcric solar Iron abundance. As a consequence of the problem in 
representing the turbulent motions on small scales in the hydrodynamical sim- 
ulations, described in ( Stefi'en et al, 2009 1, there is a smaU slope of abundance 
with equivalent width of the line, which amounts to 0.02 dex/pm. 
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Table 3. Atomic data and derived abundances for our selection of 15 Fen lines 



A 


X 


log 


9f 






yl(Fe) (HI MB) 






nm 


eV 


H 


MB 


3D (I) 


3D (F) 


(3D> (I) 


(3D> (F) 




2 84 


—2 94 


9 


7.63 


7.64 


7.60 


7.61 


7.49 


7.50 


7.56 


7.57 




z.oo 


—3 21 


— 3 21 


7.54 


7.54 


7.48 


7.48 


7.43 


7.43 


7.45 


7.45 




z.oy 


— o.oy 


fin 


7.45 


7.46 


7.40 


7.41 


7.39 


7.40 


7.39 


7.40 




3 22 


—2 23 


—2 18 


7.61 


7.56 


7.58 


7.53 


7.45 


7.40 


7.51 


7.46 


526.4812 


3.23 


-3.25 


-3.13 


7.67 


7.55 


7.68 


7.56 


7.58 


7.46 


7.64 


7.52 


541.4073 


3.22 


-3.50 


-3.58 


7.45 


7.53 


7.42 


7.50 


7.40 


7.48 


7.40 


7.48 


552.5125 


3.27 


-3.95 


-3.97 


7.45 


7.47 


7.40 


7.42 


7.42 


7.44 


7.39 


7.41 


562.7497 


3.39 


-4.10 


-4.10 


7.44 


7.44 


7.45 


7.45 


7.42 


7.42 


7.44 


7.44 


643.2680 


2.89 


-3.50 


-3.57 


7.47 


7.54 


7.47 


7.54 


7.39 


7.46 


7.44 


7.51 


651.6080 


2.89 


-3.38 


-3.31 


7.67 


7.60 


7.69 


7.62 


7.56 


7.49 


7.64 


7.57 


722.2394 


3.89 


-3.36 


-3.26 


7.59 


7.49 


7.61 


7.51 


7.55 


7.45 


7.59 


7.49 


722.4487 


3.89 


-3.28 


-3.20 


7.62 


7.54 


7.57 


7.49 


7.57 


7.49 


7.54 


7.46 


744.9335 


3.89 


-3.09 


-3.27 


7.25 


7.43 


7.29 


7.47 


7.21 


7.39 


7.27 


7.45 


751.5832 


3.90 


-3.44 


-3.39 


7.56 


7.51 


7.54 


7.49 


7.52 


7.47 


7.52 


7.47 


771.1724 


3.90 


-2.47 


-2.50 


7.54 


7.57 


7.56 


7.59 


7.43 


7.46 


7.50 


7.53 


(A(Fe)) 








7.53 


7.53 


7.52 1 7.51 


7.46 1 7.45 


7.49 1 7.48 



H MHannaford et af.U1992ll 

MB: |Melendez and Barbuy| | [2009ll 



The abundances derived from the ID models depend on the choice of the 
microturbulence parameter. When using the standard microturbulence value of 
Cmic = l.Okms"^, and the log gf values of [Melendez and Barbuy| (|2009p . the 
abundances are ^(Fe)(3D) = 7.45 ± 0.03, A(Fe)uM = 7.48 ± 0.03 for disc centre, 
and ^(Fe)(3D) = 7.48 ± 0.05, ^(Fe)HM = 7.50 ± 0.05 for integrated disc. If 
we take instead the lower microturbulence values predicted by the 3D model, 
Cmic ~ 0.7kms~^ for disc centre, and ^mic ~ 0.9kms~^ for the integrated disc 
( [Steffen et ai, 20"09| , the resulting Iron abundance is slightly higher by 0.04 
and 0.02 dex, respectively. 

From our analysis, the recommended solar Iron abundance is A(Fe) = 7.52 ± 
0.06. 

5.9. Europium 

Basing our 3D LTE analysis ( [Mucciarelli et ai, 2008[ ) on the five Full lines 
from ILawler et aP (|200ip . we obtain A(En)= 0.52 ± 0.02. The isotopic ratio 
of ^^^Eu/ (-'^^-'^Fu + -'^^'^Fu) derived from disc centre and integrated disc spectra is 
(49 ± 2.3)% and (50 ± 2.3)%, respectively. We are in perfect agreement with the 
results of lLodders et aU^UU^ . 

5.10. Hafnium 



For the abundance determination of Hafnium ( Caffau et ai, 2008bp , we consid- 



ered four Hfii lines, suggested bv ILawler et aLl(|2007p . The lines are weak and 
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blended, making the abundance determination a difBcult task. Our abundance 
determination, A(Hf)= 0.87 ± 0.04, is in close agreement with the photo- 
spheric value of ILodders et aL\ (j2009p , but larger than their meteoritic value 
of 0.73 ±0.02. 

5.11. Osmium 

Osmium and Iridium are among the heaviest stable elements and as such are 
reliable reference elements for measuring the decay of the radioactive nuclei ^ssy 
and ^■^^Th. Therefore Os is important in radioactive cosmochronology. 

The determination of the solar Os abundance is not an easy task: the suitable 
Os I lines lie mostly in the crowded violet and near UV range, and the transition 
probabilities have been revised several times in the last decades. Improved centre- 
of-gravity wavelength are given bv llvarsson et a/.1(|2003p . who recall also that Os 
has seven stable isotopes; the even isotopes which account for more than 80% 
of all Osmium have no hyperfine structure. However, according to the analysis 



of Quinet et al. (j2006|) . the isotopic splitting and the hyperfine structure of the 
odd isotopes have no influence on the Os abundance determination. 

Ilvarsson et al.\ (|2003p determined a new set of log gf for 18 Osi lines and 
compared them with previous values, finding a good agreement of their experi- 
mental results with previous determinations for the three strongest lines (290.9, 
305.8, and 330.1 nm). However, a substantial disagreement is found for other 
lines with offset and dispersion increasing for the weaker lines. The comparison 
with theoretical results shows an even larger scatter. Ilvarsson et al.\ (|2003p did 
not study the solar Os abundance, but applied their new log gf values to the 
analysis of the metal poor giant CS 31082-001, previously analysed bv lHill et al.\ 

The solar abundance of Os has been determined by [Jacoby and Aller] (11976^ 
from three Osi lines (305.8706, 330.1579, and 442.0460nm), using log gf from 
ICorliss and Bozm"anl (|1962p . to obtain A(Os)=0.70. Owing to the large errors 
often present in the ICorliss and BozmanI (|1962p data, a new determination of 
the lifetimes of nine transitions of Osmium, including the three lines used by 
[Jacoby and Aller| (|1976p . has been done by IKwiat kowski et aV (1984). These 
new log gf are systematically lower than those by Corliss and Bozman^ (I1962p 
(by about 0.5 dex) leading to A(Os)= 1.45 ± 0.10 from the equivalent widths of 
nine lines measured in the solar atlas of |Delbouille, Roland, and Neven| (|1973p . 
and using the HM model. This abundance agrees within one a with the Os 



abundance in meteorites of A(Os)= 1.37 ± 0.03 ( Lodders et al, 2009 ). 

The most recent transition probabilities and solar Osmium abundance are 
due to |Quinet et al.\ (|2006p . who also give all of the possible contaminations of 
each Os i line in the Sun or in metal-poor stars. New log gf are determined by 
these authors for 11 Osi lines and the difference with respect to the previous 
results is less than 0.1 dex for all but the two lines at 397.2 nm and 442.0 nm. 
Their newly determined log gf value for the 442.0 nm line differs by 0.33 from 
the previous value; the new log gf value of —1.20 produces an Os abundance 
in better agreement with that from other Osi lines in the star CS 31082-001. 
They determine the Os abundance using SYNTHE ( |Kurucz, 1993| ) with HM 
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Table 4. Os I lines used in the analysis with the abundance obtained from 
the EWs of disk-centre using 3D, (3D}, and IDlhd model atmospheres. The 
microturbulence was set to 1.0 km s"^ for the ID models. 



A 


EW 


X 


log 9f 




A{Os) 




Weight's 


nm 


pm 


cV 




3D 


m 


IDlhd 




305.8660 


1.03 


0.00 


-0.41<3 


1.070 


1.069 


1.020 


1 


326.7945 


0.39 


0.00 


-1.09'3 


1.211 


1.228 


1.177 


1 


330.1559 


0.72 


0.00 


-0.74<3 


1.157 


1.165 


1.114 


3 


375.2524 


0.50 


0.34 


-0.98'3 


1.442 


1.451 


1.406 


1 


397.7231 


0.07 


0.64 


-1.94'3 


1.777 


1.789 


1.744 


1 


409.1817 


0.06 


0.76 


-1.66^ 


1.533 


1.544 


1.499 


3 


442.0468 


0.23 


0.00 


-1.53^ 


1.195 


1.218 


1.165 


1 


455.0410 


0.04 


1.84 


-0.71^ 


1.415 


1.412 


1.371 


3 


463.1828 


0.01 


1.89 


-1.19^ 


1.332 


1.329 


1.288 


1 


(A(Os)> 








1.348 


1.356 


1.309 





Weight accordin g to|Quinet et al.\ II2006II 



Quinet et ar]ll2006r 
Kwiatkowski a^.1lll984ll 



and MARCS models, but do not specify whieh observational data are used. 
They discarded not only the heavily blended lines, but also the weakest ones 
(EW< 0.2 pm) because they are faint, sensitive to blends, and more affected 
by uncertainties in the log gf values. In the end, Quinet et al. (j2006p retained 
only three lines (326.7945, 330.1559, and 442.0468 nm) that give the abundances 
A(Os)= 1.40, 1.30, and 1.15, respectively, with the HM model, and A(Os)= 
1.30, 1.25, and 1.10 with a MARCS model. Their final Os abundance is A(OS)= 
1.25 ±0.11. 

We adopted the equivalent widths of the Os i lines from IKwiatkowski et aL\ 
(|1984p and derived the abundances interpolating in theoretical curves-of-growth 
(see TableE]). We adopted the partition function from lAlleiil(|1976p . If we weight 
the abundances from the individual lines as in their analysis, we obtain an 
abundance of A(0s)3d = 1-36 ± 0.19. The straight average gives yl(0s)3D = 
1.35 ± 0.21. The log gf values that we use are from [Quinet et aL\ (j2006p for all 
lines except for the 409 nm, 442 nm, and 463 nm lines for which we use the log gf 
values from IKwiatkowski et aL\ (|1984p . If we restrict our selection to the three 
lines considered most reliable by [Quinet et al. (I2006P (330.1 nm, 409.1 nm, and 
455.0 nm) we obtain a very similar value for the average (A(Os)=1.37) with an 
improved but still very large line-to-line scatter ((7=0.16). For this reason we 
recommend the result from the complete sample of lines. 

As an additional check, we fitted the profile of the line at 330.1 nm, which is 
the cleanest of our lines, with the (3D)-|-SYNTIIE synthetic spectrum, allowing 
us to take into account all of the blending lines present in the region. The abun- 
dance we find from this fit is in very good agreement with the ID result derived 
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0.0 L , : : , , : ^ : , : : ^ , : : J 

330.12 330.14 330.16 330.18 
X (nm) 

Figure 2. The 3D synthetic profile the Osl fine at 330.1 nm with yl(Os)=1.17 (green dashed 
line) is over-imposed on the disc centre solar atlas of [Delbouille, Roland, and Neven] Ill973l l 
(black solid line). 



from the equivalent width. As a further ihustration, we show in Figure[2] the 
observed profile of this line over-plotted with the 3D synthetic profile. Without 
the blending lines included, it is difficult to place the continuum in this case. 

5.12. Thorium 

The only line of Thorium in the solar spectrum that can be used for abundance 
determination is the Thii resonance line at 401.9 nm ( [Caffau et ai, 2008b[ ). This 
weak absorption lies on the red wing of a stronger Fe-Ni feature and is blended 
with weaker lines. We performed a line profile fitting analysis with a 3D syn- 
thetic profile, in order to take into account the convection-induced asymmetry, 
of the Fe-Ni blend. Giving higher weight to the results from the disc centre 
data, which resolve the line profile better, we obtained A(Tli)= 0.08 ± 0.03 
dCaffau et al., 2008b I. This result is in perfect agreement with the meteoritic 



value dLodders et al., 2009[ ). 



6. Conclusions 

The abundances of the solar photosphere, based on our CO^BOLD solar model 
atmosphere, are reported in TableO We find that 3D corrections are small in 
general, and not responsible for the systematic lowering of the solar abundances 
over the past years. With the abundances given in TableO and the abundances 
from ILodders et aD (|2009p . photospheric where available, for the elements that 
we did not analyse, we obtain a solar mctallicity oi Z = 0.0153, and Z/X = 
0.0209. This value can be compared to Z = 0.0141, and Z/X = 0.0191 of 
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Table 5. The recommended abundances of the solar photosphere for the elements considered 
in this work. 



ElGmGiit 


Ion 

state 


AbuiidSjiicc 


N 
lines 


RefGrencG 


Li 


I 


1.03 ±0.03 


1 




C 


I 


8.50 ± 0.06 


45 


ICaifau et ai. IliOlOII 


N 


I 


7.86 ±0.12 


12 


ICaffau et aLlll2009ll 


O 


I 


8.76 ±0.07 


10 


ICaffau et a;.lll2008all 


P 


I 


5.46 ± 0.04 


5 


ICaffau et al. (2007b') 


s 


I 


7.16 ±0.05 


7 


ICaffau and Ludwigi lliiOOTIl ICaffau et aZ.lllM07all 


K 


I 


5.11 ±0.09 


6 




Fc 


II 


7.52 ±0.06 


15 




Eu 


II 


0.52 ±0.03 


5 


IMucciarelh et aZ. 1(1200811 


Hf 


II 


0.87 ±0.04 


4 


ICaffau et aLI(l2008bll 


Os 


I 


1.36 ±0.19 


9 




Th 


II 


0.08 ± 0.03 


1 


ICaffau et aLI (I2008bll 



ILodders et all mm and Z = 0.0134, and Z/X = 0.0181 of 

mm . 

Our spectroscopic abundances of O and Fe are in agreement, within the 
indicated mutual errors, with what is derived from hehoseismic constraints: 
A{0)= 8.86 ± 0.045, A{Fe)= 7.50 ± 0.048 pelahaye and Pinsonneault, 2006[ ). 
There is also a good agreement for the key elements C, N and O, with the most 
likely spectroscopic abundances recommended by [Pinsonneault and Delahaye| 
A{C)= 8.44 ± 0.06, ^(N)= 7.96 ± 0.10, ^(0)= 8.75 ± 0.08. 

We conclude that our derived solar metallicity goes in the direction 
of reconciling atmospheric abundances with hehoseismic data. According to 
[Turcotte and Wimmer-Schweingruber| (|2002p . the abundance of metals in the 
solar photosphere may have decreased by up to 0.04 dex by diffusion during the 
solar lifetime. This would improve the correspondence between spectroscopic 
and hehoseismic results even further. 
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